Towards a variational principle for motivated vehicle motion 
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We deal with the problem of deriving the microscopic equations governing the individual car 
motion based on the assumptions about the strategy of driver behavior. We suppose the driver be- 
havior to be a result of a certain compromise between the will to move at a speed that is comfortable 
for him under the surrounding external conditions, comprising the physical state of the road, the 
weather conditions, etc., and the necessity to keep a safe headway distance between the cars in front 
of him. Such a strategy implies that a driver can compare the possible ways of his further motion 
and so choose the best one. To describe the driver preferences we introduce the priority functional 
whose extremals specify the driver choice. For simplicity we consider a single-lane road. In this case 
solving the corresponding equations for the extremals we find the relationship between the current 
acceleration, velocity and position of the car. As a special case we get a certain generalization of the 
optimal velocity model similar to the "intelligent driver model" proposed by Treiber and Helbing. 
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I. INTRODUCTION 

The fundamentals of the traffic flow dynamics are far 
from being well established because up to now it is not ev- 
ident what the specific form of the microscopic equations 
governing the individual car motion should be. The mat- 
ter is that the car motion is controlled by the motivated 
driver behavior rather than obeys Newton's laws. In fact, 
the behavior of a driver is due to a certain compromise 
between the will to move at a speed that is comfortable 
for him and attained on the empty road, on one side, and 
the need to avoid possible traffic accidents, on the other 
side. So, comparing the car ensembles with physical sys- 
tems it is not obvious beforehand that there is a direct 
relationship between the acceleration of a given car and 
the positions and the velocities of the other cars as is the 
case for physical particles. In addition, the linear super- 
position typical in the interaction of physical particles is 
not self-evident to hold for the vehicle interaction. 

By contrast, on macroscopic scales the car ensembles 
exhibit a wide class of critical and self-organization phe- 
nomena widely met in physical systems (for a review see 
Ref. |l| ||, ^, ^). It should be pointed out that fish and 
bird swarms, colonies of bacteria, pedestrians, etc. also 
demonstrate similar cooperative motion (for a review see 
Ref. H and also Ref. ||, H, 0)- So the cooperative behavior 
of many particle systems, including social and biologi- 



cal ones, seems to be of more general nature then the 
mechanical laws and to find out what microscopic regu- 
larities are responsible for the cooperative phenomena in 
the general case is a challenge problem. The traffic dy- 
namics is widely studied in this context also due to the 
large potential for industrial applications. 

The currently adopted approach to specifying the mi- 
croscopic governing equations of the individual car mo- 
tion is the so-called social force model or the generalized 
force model. Its detailed motivation and description can 
be found in Ref. |[ ^, here we only touch its basic 
ideas. At each moment t of time a given driver a gets up 
or lets down the speed Va of his car, or keeps its value 
unchanged depending on the road conditions and the ar- 
rangement in the neighboring cars: 
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Here the term fa(va) taken typically in the form 
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describes the driver tendency to move on the empty road 
at a certain fixed speed depending on the physical 
state of the road, weather conditions, the legal traffic 
regulations, etc. The relaxation time Tq characterizes 
the acceleration capability of the given car as well as the 
delay in the driver control over the headway. The term 
faa'{xaiVa\xa' ,Vai) dcscribes the interaction of car a 
with car a' [a' ^ a) that is due to the necessity for 
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driver a to keep a certain safe headway distance between 
the cars. The force faa'iXa,Va \xa',Va') is assumed to 
depend directly only on the velocities Va, Va' and the 
position Xa , Xa' of this car pair and being of nonphysical 
nature does not meet the third Newton's law, i.e. in 
the general case faa' ^ —fa'a- Paper \ ^a\ presents and 
discusses possible generalized ansatzes of the dependence 

faa' (*^a ; '^a \ -^a' ; '^a' ) ■ 

Special cases of this model bear own names. In par- 
ticular, for a single-lane road when all the cars can be 
ordered according to the position on the road in the car 
motion direction, 

■ ■ • < Xa-2 < Xa-l < Xa < Xa+1 < Xa+2 < ■ ■ ■ , 

solely the interaction of the nearest neighboring cars a 
and a+1 is taken into account, i.e. faa' 7^ for a' = a-\-l 
and, may be, a' = a — 1 only. For this case Bando et 
al. Ilj, ^ proposed the optimal velocity model that 
describes the individual car motion as 
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where i?opt(A) is the steady-state velocity (the optimal 
velocity) chosen by drivers for the given headway dis- 
tance A = Xa+i — Xa between the cars. This model 
and its modification successfully were used to explain the 
properties of the "stop-and-go" waves developing in dense 
traffic on single-lane roads (see, e.g. Ref. [l 4, | l 5[, p^ , [T^ , 
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However, on multilane highways the behavior of traffic 
fiow becomes sufficiently complex because of the strong 
correlations in the car motion at different lanes (for a re- 
view see Ref. |l|, mi). In this case it is not sufficient 
to confine the consideration to the interaction between 
the nearest neighboring cars only and one has to spec- 
ify several independent components of the social forces 
{faa'{xa,Va\xa' ,Va')}- As a rcsult thc number of es- 
sential fitting parameters entering the social forces in- 
creases substantially. A possible way to overcome this 
problem is to formulate a mathematical principle char- 
acterizing the strategy of driver behavior in terms of 
a certain functional quantifying the objectives pursued 
by drivers. Constructing this functional may be also 
at the phenomenological level can be lightened by the 
clear physical meaning of the driver objectives. Then, 
using standard techniques one will derive the required 
governing equations. The given work presents our first 
steps towards such an approach. It should be noted 
that the derivation of microscopic governing equations 
for systems with motivated behavior based on a cer- 
tain "optimal self-organization" principle was discussed 
in Ref. |, ^ |3^, ^ The main idea of this approach 
is the assumption that individuals try to minimize the 
interaction strength or, what is actually the same, to 
optimize their own success and to minimize the efforts 
required for this. 

Our approach is related to the concepts of mathemat- 
ical economics, namely, to the concepts of preferences 



and utility (see, e.g. Ref. |4|). We suppose that at each 
moment of time a driver plans his motion in a cer- 
tain way in order, first, to move as fast as possible and, 
second, to prevent traffic accidents. In particular, in 
the previous paper using this idea we developed a 
model explaining in a simple way the experimentally ob- 
served sequence of the first order phase transitions from 
the "free flow" to the jam phase through the "synchro- 
nized mode" |3^, |8| . Such a strategy actually implies 
that a driver evaluates any possible path of his further 
motion, {x(Oj ^ > ^0}) with respect to its preferabil- 
ity. In other words, a driver can compare any two paths 
{Xi(0) i > to} and {x2{t), t > to} and decide which 
of them is more preferable, for example, Xi(^)- The lat- 
ter relation will be designated as Xi(0 ^ X2(i)- Ob- 
viously, the given relation exhibits the transitivity, i.e. 
if Xiit) >- X2{t) and X2{t) >- Xsit) then ^ Xsit)- 

In this case we may seek for a priority functional C {x} 
meeting the condition C{xi} > C{x2} when and only 
when xi ^ X2- Finally the driver chooses the best path 
Xopt{t) of his further motion maximizing the priority 
functional C{x}- So its extremals have to satisfy the 
desired microscopic governing equation of individual car 
motion. It should be noted that chosen path Xopt{t) of 
the planned motion specifies the acceleration at the cur- 
rent time moment to rather than the real trajectory x(t) 
of the car motion because at the next time t > to the 
driver again plans his motion in the same way, introduc- 
ing the corrections caused by changes in the surround- 
ings. The same concerns the car velocity and its accel- 
eration, therefore to avoid possible misunderstanding we 
will designate the real velocity and acceleration of the 
car as v{t) and a(t) whereas these values corresponding 
to the optimal Xopt{t) will be labelled by {i^it), t > Iq} 
and {^{t), t > to}, respectively. 

In this way the problem of specifying many indepen- 
dent components of the social forces is replaced by con- 
structing the priority functional describing the driver 
compromise between the will to move as fast as being 
allowed by the physical state of the road and the neces- 
sity to avoid possible traffic accidents. So to obtain the 
priority functional we may apply to the general assump- 
tions about the driver behavior. 



II. VARIATIONAL PRINCIPLE FOR THE 
INDIVIDUAL CAR MOTION 

At the first step we should determine the collection of 
phase variables characterizing the quality of a given car 
motion. We note that for the driver under consideration 
the neighboring car arrangement and its evolution should 
be regarded as given beforehand. Indeed, it cannot be 
directly controlled by him and, so, has to be treated as 
the external conditions. Formulating this problem we 
actually assume the existence of a certain collection of 
variables taken at the current time moment t that com- 
pletely quantify the priority measure of the car motion 
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at the same time. Adopting the latter assumption we 
may construct the priority functional H, {x} in terms of 
a certain integral of a function T of the phase variables 
with respect to time. 

Applying to the conventional driver experience we will 
characterize the individual car motion at each time mo- 
ment t by its position on the road x(£), the velocity 
and the acceleration a(t\ For a multilane highway, for 
example, the position of a car should also bear the in- 
formation about the lane occupied by the car, but this 
problem will be considered elsewhere and in the present 
paper we confine our consideration to a single-lane road 
only. In this property a car is distinct from a physical 
particle because the motion of the latter is completely 
determined by the current position and velocity. The 
variables x(t\ v(t), and a{t), however, exhibit different 
behavior. The coordinate x{t) and the velocity v{t) of 
the car vary continuously, i.e. the driver cannot change 
them immediately. In contrast, the acceleration a{t) may 
exhibit sharp jumps because it is the acceleration that is 
controlled directly by the driver without remarkable de- 
lay. In such an analysis it is quite reasonable to ignore 
the short physiological delay in the driver behavior to 
changes in the surroundings, allowing sharp jumps in the 
dependence a{t). Therefore, planning his further path of 
motion the driver regards the position xq :— x{to) and 
the velocity vq :— w(to) of the car at the current moment 
of time ^0 as the initial data. 

Now let us write the general form of the priority func- 
tional £ {x} for a trial path {x(i), t > to} of the further 
motion: 

oo 

C{x} = - J dt exp(^-^^-j^^T{t,x,v,n7), (3) 
to 



where J-^(t, x, w) is the density of the path priority mea- 
sure, V := dx/dt and w := d'^x/^t'^j and the exponential 
cofactor reflects the fact that drivers can monitor the 
traffic flow state and so plan the motion only inside a 
certain region of length £ in front of them. In addition, 
the leading minus has been chosen in order to reduce 
finding the maximum of the functional C {x} to getting 
the minimum of integral (|^) as it is the typical case in 
physical theories. The direct dependence of the func- 
tion T(t,x, VtVj) on the time t reflects the effect of the 
surroundings, i.e. the physical road state and the neigh- 
boring car arrangement on the driver planning. 

According to the adopted assumption the driver 
chooses the path Xopt {t) of his further motion that maxi- 
mizes the functional C{t, Xi ^) and at the current time 
io together with other all trial paths meets the 

conditions 

x{h) = 2^0, ^{to) = vo . (4) 

Besides, the present paper analyzes the car motion in 
the traffic flow, i.e. does not consider any way by which 
a fixed car leaves the traffic fiow, for example, to stop. 
So we assume all the trial paths to exhibit bounded vari- 
ations, i.e. there are such constants , , , and 

C'jt - to) < x{t) < C'^{t-to), 

Kt)<a, \wit)\<C^. (5) 

Then using the standard technique we get the govern- 
ing equation for the extremals of the priority functional 
C {x} following from the condition 6C {x} = at x(<) = 
Xopt(t): 
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T{t,x-v.,^)^ = 0. 



(6) 



By virtue of (|^) the function x, J^, 'S?) exhibits bounded variations, which enables us to integrate equation (|^) 
twice with respect to time t reducing it to the following 
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Equation ([7|) relates the planned acceleration ru to the car 
position X and the velocity v. Subjecting this equation 
to the initial conditions (Q) we can find the optimal path 
Xopt(i) and then setting t — to we will get the desired 



relationship of the real current position x{to) and veloc- 
ity v{to) of the car with the acceleration a{to) that driver 
takes under these conditions. Exactly the expression to 
be obtained in such a way is the microscopic governing 
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equation of the individual car motion. In addition, it 
should be noted that equation (||) is of fourth order as 
it must because a trial path is fixed in part by the posi- 
tion and the velocity at the initial and terminal points. 
However, in the case under consideration the characteris- 
tics of the terminal point are replaced by conditions (|^), 
allowing us to reduce the order. 

Now let us demonstrate the proposed approach ana- 
lyzing a simple example. 



where the time scale r characterizes the acceleration ca- 
pability of the given car. The driver monitoring the car 
arrangement in front of him can predict the situation 
development, which is described in terms of the linear 
dependence of the optimal velocity t?opt(i, x) on the time 
t and the distance x: 



opt 



l + £t- 



t-tn 



Xq 



(9) 



III. THE GENERALIZED OPTIMAL VELOCITY 
MODEL 

In constructing the priority functional we assumed that 
in the steady-state traffic flow a driver prefers to move at 
a certain speed i?opt(i, x) depending on the surroundings 
and given beforehand. Besides, we consider the car mo- 
tion without acceleration to be the best way of driving. 
Therefore we adopt the following ansatz: 



T{t,x,v,a) = ^ b 



(8) 



where z^^pj = i}opt{to,XQ) and et. Ex are constants re- 
garded here as small parameters of the same order. In 
addition, the difference i'(t)/'i?opt — 1 is also assumed to 
be of order of et ^ Sx- We have adopted the linear de- 
pendence of ■!?opt(i, x) on t and x because it seems quite 
reasonable that a driver uses the linear approximation in 
estimating the position of the cars in front of him. 

Substituting (^) and (||) into ^ and truncating all the 
terms whose order exceeds e* ~ e^; ~ (^(O/^^opt ~ 1) 
get 



T'^w{t) + j dt' 



exp 



opt J 



(10) 



Multiplying equation (|o|) by the factor exp(— t^Qpt^/f), 
differentiating the obtained result with respect to t, and 
solving the equation subject to the initial condition (Q) 
we get the dependence v(t) obeying (|l0|). The differen- 
tiating v(t) at the current moment of time point we 
obtain the desired microscopic governing equation for the 
individual car motion: 



where the constants 



(11) 



and a 



opt 



2i 



In particular, let the optimal velocity 'do-pt{t,x) = 
i?opt(A) be specifled completely by the headway distance 
A — Xa+i — Xa between the given car a and the nearest 
one a -|- 1 in front of it. Then within the linear approxi- 
mation of the situation development the driver of the car 
a can anticipate that the headway distance will change 
in time as 

A(t) = A{to) + K+i(to) - vMKt - to) . 

This expression together with the dependence '!?opt(A) 
enables us to calculate the speciflc value of the constant 



Et and then to rewrite formula ( pi] ) as 



dv _ 1 



V - 'i?opt(A) - kt6v 



rf^opt(A) 
dA 



(12) 



where we have introduced the relative velocity Sv — 
Va+i — Va of the car a -I- 1 with respect to the given car 
a and omitted the argument to assuming all the values 
to be taken at the current moment of time. 

It should be noted that the obtained expression (^2|) is 
similar to the phenomenological dependence i^opt(A,'5w) 
generalizing the standard optimal velocity model , the 
"intelligent driver model" proposed by Treiber and Hel- 
bing 1^ (see also Ref. ^o|) . 



IV. CONCLUSION 

Concluding the present paper we once more remember 
its main key points. 

We deal with the problem of deriving microscopic equa- 
tions governing the motion of individual vehicles. The 
currently adopted approaches similar to the social force 
model relate in the spirit of Newton's laws the accelera- 
tion of a given car to the position and velocities of the 
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neighboring cars. In order to apply such models to anal- 
ysis of the traffic dynamics one has to specify all the es- 
sential components of the corresponding effective forces 
acting between the cars. However when the vehicle in- 
teraction becomes sufficiently complex as it is the case, 
for example, for dense traffic on multilane highways such 
an approach meets the problem of large number of fitting 
parameters. 

The present paper proposed a possible way to avoid 
the aforementioned difficulty. The main idea is to de- 
scribe at the first step the strategy of the driver behavior 
determined by the compromise between the driver will to 
move as fast as possible on the given road, on one hand, 
and the necessity to keep a safe headway distance and 
not to interfere with cars moving at neighboring lanes, 
on the other hand. This assumption actually implies 
that a driver can compare various ways of his further 
motion with respect their relative preference and choose 
the best (optimal) one at each time moment. This choice 
gives the relationship between the acceleration of the car 



under consideration and the arrangement of neighboring 
cars. 

Following the concepts of mathematical economics we 
have introduced a priority functional in order to quan- 
tify the driver choice. The extremals of this priority 
functional describe the optimal path of the driver fur- 
ther motion. At the present paper we have considered 
traffic flow on a single lane road, constructed in this case 
the general form of the priority functional, and derived 
the equations for its extremals corresponding to the car 
motion with traffic flow. The latter means that here we 
do not analyze the paths by which a car enters or leave 
traffic flow on the given road because the question de- 
serves an individual investigations. 

By way of example, we have considered a special case 
leading to an expression relating the current acceleration 
of a fixed car to the headway distance between this car 
and one in front of it as well as their relative velocity. 
The obtained equation turned out to be similar to the 
"intelligent driver model" by Treiber and Helbing. 
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